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The ratio R T)1 (T(1S)) = r r ^ ls - j ^ T + T ~ /r r ( 1S )^ M + M - is measured using a sample of (121.8 ± 
1.2) x 10 6 T(3S) events recorded by the BABAR detector. This measurement is intended as a test of 
lepton universality and as a search for a possible light pseudoscalar Higgs boson. In the standard 
model (SM) this ratio is expected to be close to 1. Any significant deviations would violate lepton 
universality and could be introduced by the coupling to a light pseudoscalar Higgs boson. The 
analysis studies the decays T(3S) — > Y(1S)-k + tt~ , T(1S) — > l + l~ , where I = h,t. The result, 
R Tt ,(T(lS)) = 1.005 ± 0.013(stat.) ± 0.022(syst.), shows no deviation from the expected SM value, 
while improving the precision with respect to previous measurements. 

PACS numbers: 13.20.Gd, 14.80.Ec, 12.60.Fr 



In the standard model (SM), the couplings of the gauge Aside from small lepton- mass effects, the expression for 
bosons to leptons are independent of the lepton flavor. the decay width T(1S) — > l + l~ should be identical for all 
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leptons, and given by 

r r (i5)^«=4a Q b ^^(l + 2^)^1-4^, (1) 

where a is the electromagnetic fine structure constant, 
Qb is the charge of the bottom quark, i?„(0) is the non- 
relativistic radial wave function of the bound bb state 
evaluated at the origin, Mr is the Y(1S) mass and Mi is 
the lepton mass. In the SM, one expects the quantity 

Ru>(T(lS)) = I rr(15) "" (2) 
1 r(is)^/'/' 

with l,V — e,/i, r and V 7^ I, to be close to one. In 
particular, the value for R Ttl (Y(\S)) is predicted to be 
- 0.992 0. 

In the next-to- minimal extension of the SM |3j , devi- 
ations of Rw from the SM expectation may arise due to 
a light CP-odd Higgs boson, A . Present data [4| do not 
exclude the existence of such a boson with a mass below 
10 GeV/c 2 . Among other hypothetical particles, A may 
mediate the following processes [l|: 

T(lS)-> A°7^/ + r 7 (3) 

or 

T(IS) -> Tfc(lS)7, -)■ A -> Z+r . (4) 

The latter implies a mixing between A and 776 (IS), 
which is a 1 S'o bb state and therefore not expected to 
decay to a lepton pair to leading order in the SM. 

If the photon is energetic enough to be detectable, 
a monochromatic peak in the photon spectrum recoil- 
ing against the lepton pair could be an indication of 
new physics (NP) 0, Alternatively, if the photon 
remains undetected, the lepton pair would be ascribed 
to the Y(1S) and the proportionality of the coupling of 
the Higgs to the lepton mass would lead to an appar- 
ent violation of lepton universality. This effect should 
be larger for decays to t + t~ pairs, and enhanced for 
higher-mass T(nS) and r)b{nS) resonances. The devi- 
ation of Rw from the expected SM value depends on 
X d = cos 6a tan (3 (where 9 a measures the coupling of 
the T(1S) to the A , and tan/3 is the ratio of the vac- 
uum expectation values of the two Higgs doublets) and 
on the mass difference between ^4° and r]b(lS). Assum- 
ing Xd = 12 (a representative value evading present lim- 
its (3), T(rj b (lS)) = 5 MeV, and M % (is) as measured 
in [7(, the deviation of R TtIl (Y(lS)) may be as large as 
~ 4%, depending on the A mass [l|. A measurement 
of this ratio has already been performed, with the result 
Rr^YilS)) = 1.02 ± 0.02(stat.) ± 0M(syst.) @]. 

This paper focuses on the measurement of R TtJ ,(T(lS)), 
in the decays Y(3S) ->■ Y(1S)tt+tt- with Y(1S) ->■ l+l~ 
and I — fi,r. In this analysis only r decays to a single 



charged particle (plus neutrinos) are considered. This 
choice simplifies the analysis; in particular it results in 
final states of exactly four detected particles for both 
the /_t + /i~ and t + t~ samples. The data collected at the 
Y{3S) resonance by the BaBAR detector at the PEP-II 
storage rings correspond to 28 fb _1 . About one tenth 
of the complete available statistics is used to validate 
the analysis method and the signal extraction procedure. 
This validation sample is discarded from the final result 
in order to avoid any possible bias. A sample of 2.4 fb _1 
collected about 30 MeV below the Y(3S) resonance (off- 
resonance sample) is also used as a background control 
sample. 

The BABAR detector is described in detail elsewhere 0, 

ED- 

The event selection is optimized using Monte Carlo 
(MC) simulated events, generated with EvtGen [rij |. 
GEANT 12] is used to reproduce interactions of parti- 
cles traversing the BABAR detector, taking into account 
the varying detector conditions and beam backgrounds. 
Final state radiation (FSR) effects are simulated using 
PHOTOS 

The selection requires exactly four charged tracks, each 
with transverse momentum 0.1 < pt < 10 GeV/c, geo- 
metrically constrained to come from the same point. The 
distance of closest approach to the interaction region of 
each track must be less than 10 cm when projected along 
the beam axis and less than 1.5 cm in the transverse 
plane. The ratio of the 2nd to 0th Fox- Wolfram moments 
(R2) [l3| is required to be less than 0.97, and the abso- 
lute value of the cosine of the polar angle of the thrust 
axis to be less than 0.96. 

A Y(1S) — > candidate is formed by selecting 

two oppositely-charged tracks, constrained to come from 
a common vertex, and it is combined with two other 
oppositely-charged tracks, assigned the pion mass, to 
construct a Y(3S) —> Y(1S)tt + it~ candidate. 

Different selection criteria are used for the Y(1S) —> 
and the Y(1S) — > t + t~ decays, because in the lat- 
ter the presence of neutrinos in the final state leads to a 
larger contamination from the background (mainly non- 
leptonic Y(1S) decays and e + e~ — > t + t~ events). The 
Y(1S) — > candidates (D^) are selected by requir- 

ing two tracks in the final state identified as muons. This 
identification is performed by using information from dif- 
ferent subdetectors, such as the energy deposited in the 
electromagnetic calorimeter, the number of hits in the 
instrumented flux return of the magnet and the number 
of interaction lengths traversed, combined in a neural- 
network algorithm. Calculated in the e + e~ center-of- 
mass (CM) frame (l6j |. the difference between the initial 
state energy and the visible final state energy is required 
to be less than 0.5 GeV, the magnitude of the dipion 
momentum (p*„.) less than 0.875 GeV/c, and the cosine 
of the angle between the two lepton candidates less than 
-0.96. For the Y(1S) t + t~ candidates (D T ), tighter 
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selection criteria are applied to reduce background. In 
these events a large fraction of the energy is not recon- 
structed, due to the presence of neutrinos; thus the dif- 
ference between the energy of the initial state and the 
energy detected in the final state, calculated in the e + e~ 
CM frame, is required to exceed 5 GeV. Further re- 
quirements are made on the magnitude of the dipion 
momentum < 0.825 GeV/c) and on the magnitude 
of the momentum of each ir (p* < 0.725 GeV/c). The 
measured difference in the energy of the T(3S) and the 
T(15) is restricted to 0.835 < AE* < 0.925 GeV. A 
boosted decision tree [TTJ is used to further reduce the 
background, based on several event shape and kinematic 
variables such as R2 and the energy of the charged tracks 
reconstructed in the events. The performance of the clas- 
sifier is assessed using MC simulations and off-resonance 
data. 

Finally, in order to select T(3S) — > Y(1S)tt + it~ candi- 
dates, the invariant mass difference AM = M(T(3S)) — 
M(T(1S)), calculated with the reconstructed tracks of 
the final state, is required to be less than 2.5 GeV/c 2 and 
the dipion invariant mass (M W7r ) to be between 0.28 and 
0.90 GeV/c 2 . 

For events with multiple candidates, the candidate 
with the value of AM closest to the nominal value Q 
is retained as the best one. It has been verified by MC 
simulations that the selection requirements do not re- 
duce the sensitivity to NP processes. Since the possible 
NP effects, with the presence of additional photons in 
the process, should be more evident in T(1S) — > t + t~ 
events, variables that are sensitive to neutral energy are 
not used in the selection. 

The final selection efficiency for the reconstructed de- 
cay chains, estimated from a sample of MC simulated 
events, are e MM = (44.57 ± 0.04)% and e TT = (16.77 ± 
0.03)% for the p + p~ and the t+t~ final states, respec- 
tively. 

An extended unbinned maximum likelihood fit, applied 
simultaneously to the two disjoint datasets D M and D T , is 



where: 



used to extract R T n = 



where AT,. 



(AW) 



indicates the number of signal events in the (D T ) 
sample. For the sample, a 2-dimensional probability 
density function (PDF) is used, based on the invariant 
dimuon mass M„+,,- and M r + C ° the invariant mass of 
the system recoiling against the pion pair, defined as: 



M^ = Js + M^-2-V^-E* 



(5) 



where ^fs is the e + e~ center-of-mass energy and E*^ 
indicates the iy + tt~ pair energy. MC simulations are used 
to check that the two variables are uncorrelated. For 
the D T sample, a 1-dimensional PDF is used, based on 
Ml%°°. (Eq.0). The likelihood is written as: 



C7 



(6) 



C % — 

'"ext 



lin, 



(7) 



fc=i 



with i = /i or r and where Ni and N,- are the sum of 
the signal and background events, observed and expected 
respectively, in each sample. Vk is the probability to 
measure a set of physical observables in the k th event, 
defined as: 



n 



N: 



^(M^-P^M^ 



bkgn ^ybkgn ( n freco \ q-jbkgfi i \ 



(8) 



and: 



€ TT N sigfj. 



e MM N. 

Nbkgr ^bkgr 



^r t ,ti{m:%™-) 



NL 



(9) 



where N^kg^ {Nbk gT ) indicates the number of background 
events in the (D T ) sample. 

The functional forms of the PDFs describing the signal 
components are modeled from the dedicated sub-sample 
consisting of one tenth of the sample. Both the 
Af r + co _ and the M,.+ n - distributions arc described by 
an analytical function approximating a Gaussian distri- 
bution function with mean value fi but different left and 
right widths, <tl,.r, plus asymmetric non-Gaussian tails 
&l,r, defined as: 



J-(x) — exp^y 



(x 



2cr L,R + ®L,r{x - A 4 ) 2 



(10) 



All the parameters (the five parameters describing the 
M„+„- distribution, along with the mean values and the 
widths of both the M re +°_ distributions) are free in the 
fit, except for oll,r in M^°_. The off- resonance sample 
is used to model the background shapes. Constants are 
chosen for the sample, and a first order polynomial 
for the D T sample, with all the parameters free in the fit. 

The result of the simultaneous fit is R Tfl — 1.006 ± 
0.013, where the quoted error is statistical only. FigureQ] 
shows the projections of the fit results for the three vari- 
ables. 

Several systematic errors cancel in the ratio, such as 
errors on the luminosity, the T(3S) production cross sec- 
tion, and the T(3S) — > Y(lS)ir + ir~ branching fractions, 
as well as systematic discrepancies between data and 
simulation in the common event selection and in track 
reconstruction efficiencies, where a possible dependence 
on the track energy has been taken into account. The 
residual systematic uncertainties are related to the dif- 
ferences between data and simulation in the efficiency of 
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FIG. 1: 1-D fit projections for M M+M - (top) and for M^ c °_ 
(middle) in the D M sample, and for M^+ c °_ (bottom) in the 
D T sample. In each plot the dashed line represents the back- 
ground shape, while the solid line is the sum of signal and 
background contributions to the fit, and the points are the 
data. 



event selection, the muon identification, and the trigger 
and background digital filters (BGF) [18(. There is also 
a systematic uncertainty on the signal and background 
yields due to the imperfect knowledge of the PDFs used 
in the fit. 

The systematic uncertainty due to the event selection 
is evaluated by comparing the shape of each variable be- 
tween data and simulation and estimating the difference 
in the efficiency. The resulting systematic uncertainty is 
1.2%. 

The systematic uncertainty related to the difference 
between data and simulation of muon identification ef- 
ficiencies, applied only to T(1S) —> events, is es- 
timated by using two independent samples: one where 
both leptons are required to be identified as muons, and 



another where exactly one final charged track is a muon. 
The ratio of the efficiencies for requiring each sample is 
determined, both on data and on simulation; the ratio of 
the two results gives an efficiency correction of 1.023 and 
a related systematic uncertainty of 1.2%. 

The systematic uncertainties due to the differences be- 
tween data and simulation in trigger and BGF efficiency 
are small both in T(1S) — > and in Y(IS) — > t + t~ 

events, and they cancel partially in the ratio. A correc- 
tion of 1.020 is needed for the e TT efficiency, together with 
a systematic uncertainty of 0.10% for T(1S) —} t + t~ 
events, while a systematic uncertainty of 0.18% is quoted 
forT(lS') — > events. The impact of the uncertainty 

in the BGF efficiency has been found to be negligible. 

The uncertainty due to the imperfect knowledge of the 
signal and background shapes used in the fit is also es- 
timated. The systematic effect from fixing oll.r in the 
signal M r + co _ PDF is estimated by varying the fixed pa- 
rameter values by ±ler and repeating the fit procedure. 
Since the correlation between the parameters is found to 
be negligible, the parameters are varied independently 
and the deviations from the nominal fit arc summed in 
quadrature, resulting in a total effect of 1.1%. The uncer- 
tainty due to the choice of the background PDF shapes 
is evaluated to be 0.22%, by using alternative parameter- 
izations. In the fit, the same M r £ co _ functional form is 
used for both the and the D T sample, ignoring the po- 
tential difference in the trigger efficiency. The systematic 
uncertainty associated with this approximation is evalu- 
ated to be 0.6%, by re-weighting the parameters for the 
M r 1 co _ distribution with the parameters obtained from 
the t + t~ data sample, and requiring the magnitude of 
the momentum of one of the final state charged tracks 
not to exceed 1 GeV/ c. 

The M%%- variable is related only to the T(35) -> 
Y(\S)tt + tt~ transition and therefore cannot distinguish 
between T(1S) —> events and other T(1S) decays or 
the Higgs- mediated events of Eq. [3] and Eq.0J While this 
ensures sensitivity to possible NP effects, T(1S) generic 
decays could be a relevant source of background in the 
D T sample because the final state is only partially recon- 
structed. The event selection heavily reduces the yield 
of the T(IS) generic decays. It is estimated using a sim- 
ulated sample of inclusive T(1S) decays, and is found 
to be - 0.4% of the Y{1S) -> t+t~ signal yield. Since 
the hadronic T(1S) decays are not well measured, the 
simulation may not be reliable and a systematic uncer- 
tainty needs to be considered. A correction factor of 
0.996, taking into account this contribution, is applied to 
the T(1S) — s- t + t~ signal yield, and a systematic uncer- 
tainty equal to 0.4% is included as well. 

The systematic uncertainty associated with the simu- 
lation of the FSR by PHOTOS is found to be negligible. 

Finally, the finite size of the MC-simulated samples 
used to determine the efficiencies gives a contribution to 
the systematic uncertainty less than 0.1% in both the 
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leptonic final states. 

The total systematic uncertainty, obtained by sum- 
ming in quadrature all the contributions, is estimated 
to be 2.2%. 

Including all the systematic corrections, the ratio R TfJl 
is found to be: 

i? TAl (T(15)) = 1.005 ± 0M3(stat.) ± 0.022(syst.). 

No significant deviation of the ratio i? TM from the SM 
expectation is observed. This result improves both the 
statistical and systematic precision with respect to the 
previous measurement @. According to jlj, and assum- 
ing values for Xd, T(rn,(lS)) and M^^g^ as previously 
stated, the present measurement excludes an A with 
mass lower than 9 GeV/ c 2 at 90% of confidence level. 
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